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Abstract 

The nonthermal phenomena in clusters of galaxies are considered in the context 
of the hierarchical model of cosmic structure formation by accretion and merging 
of the dark matter (DM) substructures. Accretion and merging processes produce 
large-scale gas shocks. The plasma shocks are expected to be collisionless. In the 
course of cluster's aggregation, the shocks, being the main gas-heating agent, gen- 
erate turbulent magnetic fields and accelerate energetic particles via collisionless 
multi-fluid plasma relaxation processes. The intracluster gas heating and entropy 
production rate by a collisionless shock may differ significantly from that in a single- 
fluid collisional shock. Simple scaling relations for postshock ion temperature and 
entropy as functions of shock velocity in strong collisionless multi-fluid shocks are 
presented. We show that the multi-fluid nature of the collisionless shocks results in 
high gas compression, reduced entropy production and modified cr v -T, M-T and Lx- 
T scalings. The scaling indexes estimated for a simple model of a strong accretion 
multi-fluid shock are generally consistent with observations. Soft X-ray and extreme 
ultraviolet photons dominate the emission of strong accretion shock precursors that 
appear as large-scale filaments. Magnetic fields, turbulence and energetic particles 
constitute the nonthermal components contributing into the pressure balance, en- 
ergy transport and emission of clusters. Nonthermal emission of energetic particles 
could be a test to constrain the cluster properties. 
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1 Introduction 



Observational cosmology boosted by recent WMAP measurements of the cos- 
mic microwave background anisotropy, high-redshift supernovae distribution, 
light element abundances and impressive progress in X-ray studies of clusters 
of galaxies is now entering the stage of the precision science (e.g. Ostriker 
& Souradeep, 2004). Clusters of galaxies being currently the largest gravita- 
tionally bounded objects are used to study cosmic structure evolution and to 
constrain the basic parameters of the concordance (ACDM) model (e.g. Bah- 
call, 2000). In the case of clusters the precision cosmology approach requires 
an account of all the physical components contributing to the energy density 
with accuracy better than 10%. Magnetic fields and energetic particles in clus- 
ters contribute to the intracluster matter (ICM) pressure at least at a level 
above a few percent and the actual contribution could be substantially higher. 
Nonthermal radiation of the clusters from radio to gamma rays and neutrino 
provide a unique way to constrain physical parameters of the clusters. 

Most of the cluster baryonic content is in the form of hot gas, likely heated by 
a hierarchy of large scale shocks. The shocks provide relaxation of accreting 
and merging matter flows in a course of the hierarchical process of cluster for- 
mation. Collisionless shock waves are the main plasma heating agent and also 
serve as an universal source of energetic charged particles and radiation. We 
shall discuss below some aspects of collisionless shocks of plasma. Gas heating 
and entropy production in collisionless shocks depends on both microscopic 
and macroscopic processes of shock relaxation. Efficient nonthermal particle 
production results in high gas compression and possibly efficient magnetic 
field generation in the shock precursor but reduces postshock ion tempera- 
ture. Physical processes in the extended shock precursor can be considered 
as a realization of gas preheating effects. That effects may be reflected in the 
observable a v -T, M-T and L x -T scalings. 



2 Collisionless shocks in clusters: magnetic fields, gas heating and 
entropy production 

Modern high resolution cluster observations in radio, IR, optical, UV, and 
X-ray bands serve as probes to study the processes in the ICM. Radio obser- 
vations of electron synchrotron emission were the first clear sign of the pres- 
ence of GeV regime relativistic electrons (positrons) in clusters. The relativis- 
ts particles are most probably related to the shock compression/acceleration 
processes. 

Energetic particles could be an essential component in clusters. Nonthermal 
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particle acceleration at shocks is expected to be an efficient process at different 
evolutional stages of clusters. Nonlinear wave-particle interactions being the 
governing process of the SNR collisionless shock formation are responsible for 
both shock heating and compression of the thermal gas as well as for creation 
of energetic particle population. 

A direct study of collisionless shock waves in a laboratory is an extremely 
difficult task. Most of the experimental data on collisionless shock physics are 
coming from space experiments. There are direct observational data on the 
shock wave structure in the interplanetary medium with clear evidences for 
ion and electron acceleration by the shocks (e.g. Tsurutani and Lin, 1985). 

Computer simulations of the full structure of collisionless shock waves describe 
the kinetics of multi-species particle flows and magneto-hydrodynamic (MHD) 
waves in the strongly-coupled system. The problem is multi-scale. It requires a 
simultaneous treatment of both "microscopic" structure of a subshock at the 
thermal particle gyroradii scale where the injection process is thought to occur, 
and an extended " macroscopic" shock precursor due to energetic particles. The 
precursor scale is typically more than 10 9 times of the microscopic scale of the 
subshock transition region. 

2. 1 Micro and macro processes in multi-fluid shocks 

In the strong enough collisionless shocks (typically of a Mach number above 
a few) resistivity cannot provide energy dissipation fast enough to create a 
standard shock transition (e.g. Kennel et al., 1985) on a microscopic scale. Ion 
instabilities are important in such shocks that are called supercritical. 

At the microscopic scale the front of a supercritical shock wave is a transition 
region occupied by magnetic field fluctuations of an amplitude 5B/B ~ 1 
and characteristic frequencies of about the ion gyro-frequency. Generation of 
the fluctuations is due to instabilities in the interpenetrating multi-flow ion 
movements. The width of the transition region of a quasi-longitudinal shock 
wave reaches a few hundreds ion inertial lengths defined as U = c/uj pi ps 
2.3 x 10 7 n~°- 5 cm. Here u P i is the ion plasma frequency and n a is the ionized 
ambient gas number density measured in cm -3 . The transition region of a 
quasi-perpendicular shock is somewhat narrower. The wave generation effects 
at the microscopic scale have been studied in some detail by hybrid code 
simulations (e.g. Quest, 1988). The large-amplitude magnetic field fluctuation 
in the shock transition region were directly measured in the interplanetary 
medium (see e.g. Kan et al., 1990). 

The reflected ions with a gyro-radius exceeding the width of the shock tran- 
sition region can be then efficiently accelerated, via the Fermi mechanism, 
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by converging plasma fluxes carrying MHD fluctuations. The efficiency of the 
upstream plasma flow energy conversion into nonthermal particles could be 
high enough providing a hard spectrum of nonthermal particles up to some 
maximal energy e+. If the efficiency of ram energy transfer to the energetic 
particles is high enough, an extended shock precursor appears due to the in- 
coming plasma flow deceleration by the fast particle pressure. The precursor 
scale L is of the order of (c/i> s h) • A* - orders of magnitude larger than the 
width of the shock transition region. Here A* is the maximal mean free path of 
a particle in the energy-containing part of the spectrum and i> s h is the shock 
velocity. We shall later refer that energetic particles as cosmic rays (CR). 

The large scale ("macroscopic") structure of a collisionless shock can be mod- 
eled by a two-fluid approach with a kinetic description of nonthermal particles 
(see e.g. Blandford and Eichler 1987; Malkov and Drury, 2001 for a review) 
or by a Monte Carlo method (e.g. Jones and Ellison, 1991). In both methods 
some suitable parameterization of particle scattering process must be postu- 
lated a priori It has been shown that the front of a strong collisionless shock 
wave consists of an extended precursor and a viscous velocity discontinuity 
(subshock) of a local Mach number that is smaller than the total Mach num- 
ber of the shock wave. The compression of matter at the subshock can be 
much lower than the total compression of the medium in the shock wave with 
allowance for high compression in the precursor. We shall refer such shocks 
later as CR-modified. 

An important predicted feature of strong shocks with efficient CR acceleration 
is a possibility to amplify an initial seed magnetic field by orders of magnitude 
(e.g. Bell & Lucek 2001; Bell 2004). CR current upstream of the strong shock 
could drive magnetic fluctuations on precursor scale length. The scale L is 
~ (c/ifeh) • A* which is expected to be above a kpc, moreover the width L > 
100 kpc for a cluster scale shock. The pre-shock scale L is 3> 10 9 times longer 
than the subshock transition region where strong magnetic field fluctuations 
are directly produced by super-alfvenic bulk plasma flows. The amplitude of 
the fluctuating magnetic field energy density Wb is of the order of the shock 
accelerated CR pressure which is in turn a substantial fraction of the shock 
ram pressure 0.5 p\ v^ h . Here p\ is the shock upstream ambient gas density. For 
typical cluster parameters that provides a pG range magnetic field amplitude 
in kilo-parsecs range scale CR-modified shock precursor. The Faraday rotation 
measure RM provided by a strong CR-dominated shock in a cluster can reach 
values of > 10 rad m~ 2 and even higher. For the case of the so-called Bohm 
diffusion model the measure RM is proportional to the maximal energy of 
the ions in the energy-containing part of the CR-spectrum accelerated by the 
shock. Faraday rotation and synchrotron-Compton emission measurements are 
used to estimate the magnetic fields in clusters (e.g. Giovannini & Feretti 2000; 
Carilly & Taylor 2002; Newman et al, 2002). 
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Electron kinetics in collisionless shocks is different from that of ions. Most 
of the observable emission comes from the electrons so they require a careful 
study. Strong shocks transfer the ram kinetic energy of the flow into large 
amplitude nonlinear magnetic fluctuations on a short scale of the transition 
region. The thermal electron velocities in the ambient medium are higher than 



the shock speed if the shock Mach number A4 S < ym p /m e , allowing for nearly- 
isotropic angular distribution of the electrons. Non-resonant interactions of 
these electrons with large-amplitude turbulent fluctuations in the shock tran- 
sition region could result in collisionless heating and pre-acceleration of the 
electrons (Bykov and Uvarov, 1999). The presence of large-amplitude waves 
in the shock transition region erodes many of the differences between quasi- 
parallel and transverse shocks, providing the electron injection mechanism to 
be similar for these shocks. 

An exact modeling of a collisionless shock structure with an account of the non- 
thermal particle acceleration effect requires a nonperturbative self-consistent 
description of a multi-component and multi-scale system including the strong 
MHD-turbulence dynamics. That modeling is not feasible at the moment. In- 
stead we will use below some simplified description of a multi-fluid strong 
shock structure with an appropriate parameterization of the extended pre- 
shock and the gas subshock. Then some predicted observable characteristics 
of the shocks can be confronted against observational data. We will discuss 
first the effects of plasma heating by modified shocks and then some specific 
predictions for possible observational tests. 



2.2 Plasma ion heating by multi- fluid collisionless shocks in clusters 

The collisionless shock relaxation processes are multi-fluid. For strong colli- 
sionless shocks in a magnetized plasma the nonthermal particle acceleration 
effect is expected to be efficient and a significant fraction of the ram pressure 
could be transferred into high energy particles (ions only for the nonrelativistic 
shocks). A strong shock of a total Mach number Ai s ^> 1 consists of a broad 
precursor of a scale of L ~ (c/v B h) • A* due to the upstream gas deceleration 
by the nonthermal particle pressure gradient and a viscous gas subshock of 
the modest Mach number .M su b ~ 3 (see for a review Blandford and Eichler, 
1987; Jones and Ellison, 1991; Malkov and Drury, 2001). 

The total compression ratio R t of a strong MHD shock modified by an efficient 
nonthermal particle acceleration can be estimated as 




7 + 1 




7- v /l + 2( 7 2 -l)Qesc/pi^ 3 h ' 
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assuming that the energy density in the shock upstream is dominated by the 
ram pressure and that the CR escape is through the cut-off momentum regime 
(e.g. Malkov and Drury, 2001). Here Q csc is the energy flux carried away by 
escaping nonthermal particles and 7 is the effective adiabatic exponent. 

The distribution function of nonthermal particles and the bulk flow profile 
in the shock upstream region are sensitive to the total compression ratio R t . 
Thus, the exact calculation of the escape flux Q csc can be performed only in 
fully nonlinear kinetic simulations. Nevertheless, some approximate iterative 
approach can be used to make the distribution function consistent with the 
shock compression. The subshock is the standard gas viscous shock of a Mach 
number -M su b- For that simplified steady-state model of a strong CR-modified 
shock the downstream ion temperature TJ can be estimated (Bykov, 2004) 
for the shock of a given velocity if R t and -M su b are known: 

T/ 2) - 0(M sub ) • where 0(_M sub ) = ~ ( 2 7g ~ c . (2) 

7g#t>sh) (7g " l)M S ub + 2 



Single fluid strong shock heating represents the limit -M su b = M s ^ 1) since 
there is no precursor in that case. In a single fluid strong shock with M s ^> 1 
and M a 3> 1 we have 

T, (2) - 2 ■ i^ii ^ sh = 1.38 ■ 10 7 v% (K), (3) 



for any magnetic field inclination (e.g. Draine and McKee, 1993). The 7 g is the 
gas adiabatic exponent, the mass per particle /i was assumed to be [1.4/2.3] • 
itlh and v s s is the shock velocity in 10 8 cm s _1 . In the single-fluid system the 
compression ratio R t = R sn b — * (7g + l)/(7g — 1) does not depend on the shock 
velocity and the Eq.(2) reduces to Eq.(3). However in the multi-fluid shocks 
the total compression ratio depends on the shock velocity and could be higher 
than that in the single-fluid case. That implies somewhat lower postshock ion 
temperatures for the strong multi-fluid shock of the same velocity and could 
be tested observationally (e.g. Bykov, 2004). 

Consider multi-fluid CR-modified strong shocks. The gas heating mechanism 
in a CR-modified shock precursor must be specified in a simplified two-fluid 
approach, providing a connection between the total compression R t and the 
gas subshock Mach number -M su b- It is worth to recall the plasma parameter 
P = M 2 a / Ml is also the ratio of the thermal and magnetic pressures. 

The gas heating mechanism in the modified shock precursor is still under 
debate (e.g. Malkov and Drury, 2001) and is usually postulated in shock mod- 
eling. There are extensive simulations of the models with adiabatic heating 
of the shock precursor (e.g. Kang et al., 2002). The Alfven waves dissipation 
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mechanism was suggested by McKenzie and Volk (1982). The turbulent gas 
heating due to acoustic instability wave dissipation could be more efficient, 
however, if the acoustic instability develops (e.g. Malkov and Drury, 2001). 

The subshock has .M S ub ^3 in the most of the simulated strong shocks with 
adiabatic heating of the precursor. The result is not too sensitive to the ion 
injection rate if the rate exceeds some threshold to provide the shock modifi- 
cation by the accelerated ions. The scaling Ai su h ~ 2.9 ■ A4®' 13 was obtained 
by Kang et al., (2002). They also found that R t oc M 3 J 4 is a reasonably 
good approximation for the two-fluid time-dependent simulations of a CR- 
modified shock with adiabatic precursor. However, if precursor heating due to 
Alfven waves dissipation or the acoustic wave instability is efficient then for 
the plasma compression in the shock precursor could become less efficient. 

Notice that the approximations for the shock compression and heating given 
above were obtained within a simplified two-fluid steady-state approach where 
the Bohm-like diffusion model was postulated. The possibility of a shock acous- 
tic instability development is still an open question. It requires a dedicated 
study of the connection between the diffusion coefficient and the local flow 
parameters. Recent detailed numerical studies of the unsteady CR-modified 
shocks by Kang et al., (2002) a priori assumed the spatial dependence of the 
diffusion coefficient k(p,x) oc n" 1 (x) to prevent the acoustic instability. Gas 
heating due to acoustic instability waves or Alfven waves dissipation would 
change their results. 

A modest heating rate is expected in the case of adiabatic plasma compression 
in the precursor. That could be the case for very high Alfven Mach number 
shocks 1 < M 2 S < M & (i.e. 1 < M s < y/P) where the Alfven wave dissi- 
pation is not efficient and the acoustic wave instability does not develop. The 
total compression ratio is R t ~ Ai 3 ^ 4 (e.g. Blandford and Eichler, 1987). Then 
the downstream ion temperature for the case of a weak upstream magnetic 
field (and Ai su b ~ 3) can be estimated from: 

T ( 2 ) i T (i) w 2 2 . v /^7 s ( 1 + f e .y\ ( 4 ) 

Here f ei = T e /T- V The latter case of M s <C y//3 would require the plasma 
parameter f3 > 100 in the shock far upstream. The scaling can be also relevant 
to the internal shocks of Ai s < 10. Following cosmological shocks simulations 
of Ryu et al., (2003) the internal shocks are expected to be rather numerous 
in the ICM. Note here that the same temperature scaling T^/T^ oc VM S 
is appropriate for strong radiative shocks (see e.g. Bouquet et al., 2000). 

Under an assumption that the main heating mechanism of the gas in the 
precursor region is due to Alfven waves dissipation, an asymptotical estima- 
tion for the total compression of the shock R t « 1.5M 3 J 8 was presented by 
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Berezhko and Ellison (1999). The estimation was obtained assuming a con- 
stant Alfven speed in the shock upstream and it is valid if M 2 ^> ^> 1 
in the far upstream flow. One may then obtain the ion temperature just be- 
hind the strong modified shock with the preshock heating due to Alfven wave 
dissipation (Bykov, 2002) from Eq.(2): 

Tf » 0.32 ■ 0(M sub ) ■ v 5 s tn^ 3/8 B%\l + f ei )~\ (5) 

The ion temperature is measured in 10 6 K. Eq.(5) is valid under the conditions 
of M s > 1 and M s > \f]5 in the far upstream flow. We have 0.32 • 0(3) wl.2 
for 7 g = 5/3. 

An important distinctive feature of the strong CR-modified shock with gas 
heating by Alfven wave dissipation (given by Eq.(5)) is the dependence of 
the postshock temperature on the number density and magnetic field of the 
incoming plasma oc (B/y/n) 3 ^. There is a remarkable difference between the 
single-fluid shock heating where oc v 2 h for M 2 S ^> 1 and that for the 
CR-dominated multi-fluid strong shocks (given by Eq.(5) and Eq.(4)) where 
T( 2 )oc< h . 

It is convenient to introduce the scaling Rt(v 8 t) oc to describe the different 
cases of strong shock heating. Then t/ 2 -* oc 0(-M S ub) • A ' ) and the index 
a = 2(1 — A) + 5. Here the index 5 approximates the velocity dependence 
of (f>(M su b) oc ff h . The subshock Mach number M su b depends, in general, on 
M s and M^. For the Alfven wave heating case (see Eq.(5)) we have A 3/8, 
the index 5 is typically < 0.25, and thus a < 1.5. In the case of high enough 
efficiency of the turbulent heating of gas due to the acoustic instability of 
a strong shock the index a may exceed 1.5 being closer to the single-fluid 
limit. The shock velocity dependence of R t (v s h) and T^ivg^) discussed above 
should result in cluster thermodynamic parameter scalings o~ v -T and other to 
be confronted against observational data (see e.g. Rosati, Borgani & Norman 
2002 for a review). The scaling relations depend on electron-ion relaxation and 
we will discuss them below. 

2. 3 ICM entropy production by multifluid accretion shocks 

Cold gas falling into the DM dominated gravitational well passes through a 
strong accretion shock. The shock is a source of gas entropy production in 
the ICM (e.g. Knight & Ponman, 1997; Tozzi & Norman, 2001; Voit et al., 
2003). The post-shock entropy K = K$T / p 2 ! 3 used in the ICM analysis and 
simulations (e.g. Bialek et al., 2001) is related to the standard thermodynamic 
entropy s through K oc exp(s/c„). In the standard scenario with a single-fluid 
accretion shock the post-shock entropy K si oc v 2 h - p[ 2 ^ (e.g. Voit et al., 2003). 
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The multi-fluid nature of the collisionless accretion shock modifies the stan- 
dard scaling relation to be 

K m{ oc vi • [it^sh)]- (1+7g) • 0(-Msub) • pf" 7g) (6) 

The compression ratio in CR-shocks is higher than in a strong single-fluid 
shock of the same velocity resulting in reduced post-shock entropy production. 
For example, in the case of Alfven heating the post-shock entropy of a multi- 
fluid shock reduces as K mi /K si ~ (15 /M&) for M a > 15 and M 2 > M a . Here 
and below in numerical estimations we assume 7 g = 5/3 though non-thermal 
baryonic component could reduce the index j g . 

Since Rt(v 8 h) and 0(-M su b) are shock velocity dependent the simple scaling 
K oc v g h • p 1 2 ^ 3 is not valid. In CR-modified shocks K m f oc v^ h ■ p[ 7s ^ or 
K m{ oc T^ a , where v = 2-(l + 7g ) • A + 5. For the case of Alfven wave heating 
the index v < 1.25 and K m f oc T - 8 assuming 7 g = 5/3. Recently Piffaretti et 
al., (2004) found that the dispersion in the observed cluster entropy profiles 
is smaller if an empirical relation K oc T 65 is used instead of the standard 
K oc T. 

Consider a simple model of smooth accretion of cold gas through a strong 
accretion shock by Voit et al., (2003). The gas of velocity v ac accretes at a rate 
M g through the shock at a radius r ac where 

M g = 4nrl cPl v ac , v 2 ac = 2GM^-\ i = 1 - rjr ta . (7) 

Here M(t) is the cluster mass and r ta is the matter turnaround radius. Then 
the entropy K mf just behind the multi-fluid shock is expressed through T> (v ac ) 
and p 2 = Rt(vac) ■ p\. In the Alfven wave heating case K mi (t) oc (MtY 1+5 ^ 3 , 
instead of K sf (t) oc (Mt) 2 / 3 in the single-fluid regime. Multi-fluid shock results 
in a slower post-shock entropy production. As we have noted above the regime 
of CR-shock compression depends on the plasma parameter (3 in the infalling 
gas. The plasma parameter (5 is currently poorly known because the interclus- 
ter magnetic fields are not well constrained. The effects of shock modifications 
are important for both the models of smooth accretion of cold gas and for 
accretion of hierarchical structures as well. 

Preheating of accreting gas by different physical processes (e.g. due to early 
star formation in protocluster region) was suggested by Evrard & Henry 
(1991), as a possible reason for the breaking of the scaling relations for pure 
gravitational cluster compression (Kaiser, 1986). The observed high metallic- 
ity of clusters at different redshifts indicates that strong starburst activity was 
highly likely at some stage. The preheating produces some initial level of gas 
entropy ("entropy floor") (see e.g. extensive simulations by Bialek et al., 2001; 
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Borgani et al., 2001). Multi-fluid strong shocks provide another natural way 
of preheating of accreting gas. 

We have discussed above the ion temperature and entropy production in multi- 
fluid shocks. However, cluster radiation spectra (X-ray and multi-wavelength) 
depend in large on the electron temperature. 

2.4 Electron temperature in the multifluid shocks 

The electron and ion temperatures eventually equilibrate in a thick postshock 
layer due to Coulomb collisions. The complete e — i Coulomb equilibration 
requires the system age t > 10 10 T^ 2 /n, where the postshock density n is 
measured in cm -3 , the ion temperature T 6 is in 10 6 K, and t is in seconds (see 
e.g. Mewe 1990 for a review). The Coulomb equilibration rate is fast enough 
to guarantee T e Tj on the cluster scale for a Hubble time. 

The complete e — i equilibration is not always the case for a local vicinity of a 
collisionless shock, however. For the internal shocks in the ICM the electron- 
ion temperature relaxation layer can be of the order of 100-300 kpc and it 
can be even larger for the external shocks. It is highly likely that the e — i 
Coulomb equilibration will be complete in the extended precursors of multi- 
fluid accretion shocks. Alfven waves and magnetic fluctuation pre-heat the 
ions to temperatures of about 0.1-0.5 keV in the extended shock precursor. 
The column depth in the precursor could be about > 10 19 cm -2 providing 
substantial e — i equilibration in the gas entering the subshock of Ai sn b ~ 3. 
Then ions will be heated in the collisionless subshock by a factor of ~ 3.5. 

The gas density and ion temperature jumps do not necessarily coincide with 
that of the electron temperature. X-ray telescopes can resolve in principle the 
extended structure with gradual electron temperature profile behind shock 
front. On the other hand in case of merger events it is more easy to identify the 
"cold fronts" - sharp contact discontinuities between cool gas bodies moving 
through a hot ICM. Such structures were discovered in A2142 (Markevitch et 
al., 2000), A3667 (Vikhlinin et al., 2001) and other clusters. 

It should be noted, however, that while the Coulomb collisions are unavoidable 
they are not the only possible electron heating process in shocks. Electromag- 
netic field fluctuations providing the ion flow relaxation in collisionless shock 
may heat electrons on rather a short scale. 

The initial electron temperature just in the shock downstream depends on col- 
lisionless electron heating (see e.g. Cargill and Papadopulos, 1988; Bykov and 
Uvarov, 1999). Non-resonant interactions of the electrons with strong nonlin- 
ear fluctuations generated by kinetic instabilities of the ions in the transition 
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region inside the shock front may play the main role in the heating and pre- 
acceleration of the electrons, as it was shown in the model by Bykov and 
Uvarov (1999). They calculated the electron energy spectrum in the vicinity 
of shock waves and showed that the heating and pre-acceleration of the elec- 
trons occur on a scale of the order of several hundred ion inertial lengths U 
in the vicinity of a viscous discontinuity. Although the electron distribution 
function is significantly non-equilibrium near the shock front, its low energy 
part can be approximated by a Maxwellian distribution. The effective electron 
temperature just behind the front, obtained in this manner, increases with the 
shock wave velocity as T c oc v b sh with b < 2. They also showed that if the elec- 
tron transport in the shock transition region is due to turbulent advection by 
strong vortex fluctuations on the ion inertial length scale (Zj), then the nonres- 
onant electron heating is rather slow (i.e. b < 0.5), but the electrons are still 
injecting into the diffusive Fermi acceleration. A highly developed vortex-type 
turbulence is expected to exist in the transition regions of strong shocks. That 
would imply that the initial T e /T; oc ^ just behind the transition region 
would decrease with the shock velocity for Ai s ^> 1. So, the electron heating is 
more efficient in the shocks of lower velocities. There are some observational 
indications that T c /Tj oc for collisionless shocks in supernova remnants 
(e.g. Raymond, 2001). The scaling is generally consistent with the interplane- 
tary shock data compiled by Schwartz et al., (1988) (see also Rakowski et al., 
2003). The degree of electron-ion equilibration in collisionless shock is found 
to be a declining function of shock speed. In the case of strong vortex-type 
turbulence in the shock transition region (a — b) ~ 1 for the Alfven heating 
case described above. 



3 Effect of the non-thermal components on a v -T, M-T, L x -T rela- 
tions 

X-ray cluster studies accompanied by multi- wavelength observations have pro- 
vided a possibility to investigate statistical and evolutional properties of ICM 
in clusters (see e.g. Forman & Jones, 1982; Sarazin, 1988; Rosati, Borgani & 
Norman, 2002 for a review). In particular, relationships between the optical 
measured galaxy velocity dispersion a v , the X-ray bolo metric luminosity Lx, 
the X-ray temperature Tx and cluster gravitating mass M were determined 
observationally. Simulated scalings are to be confronted against observational 
data (see e.g. Arnaud & Evrard 1999, Mohr, Mathiesen & Evrard, 1999, Wu 
et al., 1999, and Rosati, Borgani & Norman, 2002 for a review) to constrain 
the models of cosmic structure evolution and to estimate the baryon fraction. 

The Lx-T and other relations were simulated by Cavaliere, Menci & Tozzi 
(1998, 1999) using the standard single-fluid Rankine-Hugoniot law. They con- 
sidered hierarchical merging history using the generalized Press-Schechter DM 
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halo distribution. That allows the authors to make predictions for the expected 
cosmological evolution of the L^-T correlation. We will discuss here some sim- 
plified models to study the effect of multi-fluid shocks on the thermodynamic 
relations. The shock velocity dependences of the compression ratio R t (v sii ) 
and post-shock temperature t/ 2 ^(i> sIi ) discussed above result in modified clus- 
ter thermodynamic parameter scalings a v -T, M-T, L x -T. 



In a simplified model of smooth accretion of cold gas by Voit et al. (2003), 
a strong accretion shock of a velocity given by Eq.(7) will have i> s h oc a v . 
Lumpy accretion of DM subhalos produce internal (merger) shocks inside the 
ICM gas heated by the external (accretion) shock. The properties of shocks 
in the large-scale structure of the Universe were simulated in the context of 
the ACDM-cosmology using PM/Eulerian hydrodynamic codes (e.g. Miniati 
et al, 2000; Ryu et al, 2003). Internal shocks with 2 < M s < 4 were found 
to be most important in energy dissipation providing ICM heating. The code 
deals with N-body CDM and single-fluid gas dynamics. However, if a strong 
accretion shock is multi-fluid, providing reduced post-shock ion temperature 
and entropy, then the internal shocks could have systematically higher Mach 
numbers 4 < M s < 8. 



Consider as a generic example accretion of a cold gas of n a ~ 10 -5 cm~ 3 and 
B ~ 0.03 — 0.1 fiG through an accretion shock located at r ac ~ 2 Mpc around 
a rich cluster of M ~ 1O 15 M . The gas velocity in the shock rest frame is 
i> s h ~ 2,000 km s _1 . Then the gas compression ratio by the strong multi- 
fluid shock R t ~ 8 and the post-shock ion temperature Tj ~ keV can be 
obtained from Eq.(5). Magnetic field can be amplified in the shock precursor 
and compressed by a sub-shock providing B ~ fiG in the post-shock flow. 
The extended multi-fluid shock precursor of a scale of a few Mpcs and of a 
width of some hundred kpcs will have a temperature of about 0.2-0.3 keV and 
magnetic field ~ 0.5 jiG. 



Further heating and compression of the accreted gas in the inner cluster will 
be provided by internal shocks. Though an internal shock of M s ~ 6 is not 
extremely strong, it may still be multi-fluid. For such a shock Eq.(5) is valid 
under the conditions of A4 S ^> 1 and M.1 ^> M. & in the far upstream flow of 
the inner shock. The later condition can be expressed as i> s g ^> 0.6 n~i\ T 7 Bz\- 
The ion temperature T 7 is measured in 10 7 K. Then considering the shocks as 
the main heating agents for the inner ICM one can estimate thermodynamic 
scaling relations. 
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3.1 a v -T scaling 

For a cluster near the virial equilibrium the internal shock velocity i> s h oc a v . 
The multi-fluid nature of CR-modified shocks would result in a v oc T][ a for 
the cluster with relaxed e — % temperatures. In the case of the shock with 
Alfven wave dissipation heating a < 1.5, and the a v -T scaling index is > 2/3. 
The statistical analysis of a sample of about 100 observed clusters given by 
Wu et al. (1999) provided the scaling a v = 10 2 - 49±a03 T a64±a02 consistent with 
a v oc Tj /3 . 



3.2 M-T scaling 

For a virialized cluster of a total mass M we have a 2 oc M 2//3 . Therefore, in 
the model of a multi-fluid shock M oc T 3//a is expected. For the case of relaxed 
e — i temperatures (on the cluster scale) that results in Mi CM oc T^ a ■ fb- Here 
the ICM gas mass Mi CM = fb • M, and fb is the baryon gas mass fraction. 
The index of Mi CM -T scaling derived by Mohr et al. (1999) from observations 
is 1.98 ± 0.18, consistent with the value > 2 predicted by the Alfven wave 
heating multi-fluid shock model (a < 1.5). The M-T relation for a sample of 
galaxy clusters, groups and elliptical galaxies was studied by Sanderson et al. 
(2003). The authors reported the M-T index 1.84 ± 0.06 for the sample. To 
apply our simple model to a sample of the objects of different cosmological 
ages one has to account for the evolution of the mean density of the Universe 
with redshift and for galaxy feedback processes (see e.g. Borgani et al., 2004). 



3.3 Lx-T scaling 

The model L x -T scaling can be estimated from the relation L x oc a\ ■ T 1 / 2 ■ 
fb ' r Xc gi ven by Wu et al. (1999). Here r Xc is the cluster core radius in 
spherical isothermal /3-model of ICM gas distribution (of index 2/3). Thus, in 
the multi-fluid shock model we have Lx oc X 4//a+ V 2 . For the preshock heating 
by Alfven waves a < 3/2 and the L x -T index is > 3.16. We did not account 
here for the weak temperature dependence of the thermal bremsstrahlung 
Gaunt-factor. Arnaud & Evrard (1999) obtained the L x -T index 2.88 ± 0.15 
from a sample of 26 clusters, while Markevitch (1998) obtained 2.64 ± 0.27 
removing photons from central regions of the clusters. It should be noted that 
contrary to the cr v -T and Micm~T relations the Lx-T relation is more sensitive 
to the assumed spherical symmetry and matter distribution inside the cluster. 
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Summary 



The multi-fluid processes of shock heating discussed above are a natural re- 
alization of the idea of accreting gas preheating. Collisionless shocks generate 
energetic non-thermal particles that can penetrate far into the upstream flow. 
The particles decelerate the flow and preheat the gas. They can also efficiently 
generate fluctuating magnetic fields of ~ /xG range on some kilo-parsecs scale. 

An important distinctive feature of multi-fluid shocks is their high gas com- 
pression Rt(v s h) that could be well above the single fluid shock limit (7 g + 
1) / (7 g — 1). At the same time entropy production for a strong multi-fluid shock 
scales as Rt(v s h)~^ s+1 ^ and it is significantly reduced compared to the single- 
fluid shock of the same velocity. The effects are due to energetic particle accel- 
eration and magnetic field generation. Energetic particles can penetrate into 
the shock upstream region. They are coupled with the upstream gas through 
fluctuating magnetic fields (including the Alfven waves). Magnetic field dissi- 
pation provides gas preheating and entropy production in the shock precursor. 
Energetic particles evacuate gas energy and momentum resulting in high gas 
compression and reduced temperature. Such a heated pre-shock region would 
appear as an extended filament of a width L ~ (c/v s h) • A* > 3 x 10 14 • ■ Bz\ 
cm. Here e* (in GeV) is the highest energy of the hard branch of the acceler- 
ated particle spectrum. Since B^q < 1 in the cluster outskirts and if the hard 
spectrum of energetic nuclei extends till ~ 10 9 GeV (c.f. Norman, Melrose & 
Achterberg, 1995) we have L > 100 kpc and even wider. A warm gas (~ 0.2 
keV) emission filament found with XMM-Newton in the outskirts of Coma 
cluster by Finoguenov et al. (2003) could be an extended heated precursor of 
a strong multi-fluid accretion shock. 

Energetic nuclei can be stored in cluster magnetic fields for Hubble times (e.g. 
Berezinsky, Blasi & Ptuskin, 1997). Nonthermal particles and magnetic fields 
contribute to the total pressure of the ICM. Nonthermal electromagnetic emis- 
sion and neutrinos due to energetic particles interactions potentially provide 
tests to constrain the ICM nonthermal components. Diffuse synchrotron radio 
halos were reliably detected in some clusters (e.g. Giovannini & Feretti, 2000; 
Carilli & Taylor, 2002) indicating the presence of highly relativistic electrons. 
There are some evidences for a presence of excesses above the dominating 
thermal emission of the ICM hot gas in both hard X-rays and EUV. Brunetti 
et al. (2004) argued that present observations of non-thermal radiation from 
clusters of galaxies can be explained even if the pressure of relativistic hadrons 
in the intracluster medium does not exceed 5-10% of that of the thermal gas. 
A reliable detection or a meaningful constrain of the hard X-ray tails with 
INTEGRAL or ASTRO-E2 and future gamma-ray observations with GLAST 
would be very important to identify the nature of the nonthermal processes 
in clusters. The existing upper limits for emission above 100 MeV obtained 
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with CGRO EGRET (e.g. Reimer et al., 2003) can not yet rule out substantial 
energy density of energetic nuclei (see also Blasi, 1999). A population of ener- 
getic nuclei of energy < GeV is hard to detect with current telescopes even 
if they have a substantial pressure comparable with that of the hot thermal 
gas. Ultra high energy CR acceleration above 10 9 GeV by cosmic shocks was 
suggested by Norman, Melrose & Achterberg (1995) as a plausible scenario. 

Study of the cooling flows spectra with RGS and spatially resolved X-ray spec- 
troscopy with EPIC aboard XMM-Newton (e.g. Kaastra et al., 2004) revealed 
significant systematic deficit of the emission of the low temperature gas as com- 
pared to the isobaric cooling flow model predictions. Magnetic fields should be 
progressively more important in the cluster central region. While the strong 
multi-fluid shocks are expected to dominate at the outskirts of the cluster the 
hierarchical clustering model predicts moderate and weak shocks in the inner 
parts of the cluster as it was demonstrated by Ryu et al., (2003). One may 
also consider a magnetic field reconnection induced by cluster substructure 
motions. In the scenario the local energy release rate due to magnetic fields 
reconnection would be oc B c oc rf a with c > 2. Dissipation of the weak shocks 
from the reconnection events and that produced by fast moving substructures 
would preferentially produce energetic particles thus allowing for the non-local 
heating effect different from the standard heat-conduction (Bykov, 2002). The 
cold-phase filling factor could be strongly reduced in that case. 
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